Determination of the Fermi Surface of MgB 2 by the de Haas-van Alphen effect 



CO 

o 
o 

CN 
< 

OO 



S3 
O 
O 

Oh' 



^3 

o 
o 



> 
in 
co 
^1- 

o 

CO 

o 

-I— > 



i 

o 
o 



X 



A. Carrington a , P.J. Meeson a , J.R. Cooper b , L. Balicas c , N.E. Hussey a , E.A. 

Yelland fc , S. Lee d , A. Yamamoto d , S. Tajima d , S.M. Kazakov 6 and J. Karpinski e 
a H. H. Wills Physics Laboratory, University of Bristol, Tyndall Avenue, BS8 1TL, United Kingdom. 
b Interdisciplinary Research Centre in Superconductivity and Department of Physics, 
University of Cambridge, Madingley Road, Cambridge CB3 OHE, United Kingdom. 
c National High Magnetic Field Laboratory, Florida State University, Tallahassee, Florida 32306, U.S.A. 
1 Superconductivity Research Laboratory, International Superconducting Technology Center, Tokyo 135-0062, Japan, and 
e Laboratorium fur Festkorperphysik, ETH Zurich, CH-8093 Zurich, Switzerland. 

(Dated: February 2, 2008) 

We report measurements of the de Haas- van Alphen effect for single crystals of MgE>2 , in magnetic 
fields up to 32 Tesla. In contrast to our earlier work, dHvA orbits from all four sheets of the Fermi 
surface were detected. Our results are in good overall agreement with calculations of the electronic 
structure and the electron-phonon mass enhancements of the various orbits, but there are some 
small quantitative discrepancies. In particular, systematic differences in the relative volumes of the 
Fermi surface sheets and the magnitudes of the electron-phonon coupling constants could be large 
enough to affect detailed calculations of T c and other superconducting properties. 



There has been rapid progress in understanding the 
physical properties of magnesium diboride since the dis- 
covery of superconductivity at 39 K just over two years 
ago. The accepted consensus is that it is an s-wave, 
phonon-mediated superconductor but with some highly 
unusual properties. The most important of these are the 
anomalously high T c and the existence of two, almost dis- 
tinct, superconducting gaps. This understanding is based 
on many experiments and on theoretical calculations of 
the unusual electronic structure of MgB 2 0, H B • 

There have been two direct experimental probes of the 
Fermi surface (FS) structure of MgB 2 ; angle-resolved 
photoemission spectroscopy || and the de Haas-van 
Alphen effect (dHvA)[f|. Our previous dHvA study 
showed good agreement between theory and exper- 
iment regarding the areas of the orbits that were ob- 
served and their electron-phonon mass enhancements. 
However, only three out of nine predicted dHvA orbits 
were observed (labelled 1,2 and 3 in Fig. ^| and so in- 
formation was only obtained about two of the four FS 
sheets predicted by band calculations. It was not en- 
tirely clear whether the unobserved orbits were missing 
simply because of their relatively short mean- free-paths, 
or whether the topology of the other sheets was substan- 
tially different. This issue clearly affects calculations of 
many physical properties (especially T c ). In this Letter, 
we report new measurements of the dHvA effect in mag- 
netic fields up to 32 Tesla for two single crystals. These 
new results give evidence for orbits on all four sheets of 
the Fermi surface. 

Experiments were conducted on two different crystals 
grown by different groups. Sample B is one of the crys- 
tals studied previously [6j and was grown in Tokyo by 
high pressure synthesis using natural boron. Sample K 
was grown in Zurich by a similar technique but isotopi- 
cally pure boron-10 was used as a starting ingredient [7J. 
Other crystals from the same batches were found to have 




FIG. 1: Calculated Fermi surface topology, with possible 
dHvA extremal orbits (for frequencies <10 kT) indicated. 
Figure adapted from Kortus et al. Q]. 



T c (onset) of 38.0 K and 37.7 K for the Tokyo and Zurich 
crystals respectively. dHvA oscillations were observed by 
measuring the torque (r) with a sensitive piezo-resistive 
cantilever technique S}. Both samples were studied ex- 
tensively at Bristol in fields up to 20.5 T and temper- 
atures down to 0.3 K, and in fields up to 32 T at the 
NHMFL. 

The first harmonic of the oscillatory part of the torque 
for a 3D Fermi liquid is given by 0, [3 
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where F is the dHvA frequency [F — (h/2ire)A, A 
is the extremal orbit area in k-space]; C is the cur- 
vature factor, 7 is the phase; Rrj, Rt and R$ are 
the damping factors from impurity scattering, temper- 
ature and spin splitting respectively. The Dingle factor, 
Rd = e xp( — ^j-), where ms is the unenhanced band 
mass 0, and r is the scattering time. An equivalent 
expression for Rd is Rd — exp(— ), which shows 
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FIG. 2: Fourier transforms of the raw high field data for crys- 
tal K at 1.4 K and fields from 20— > 32 T. (the square root 
of the FFT amplitude is plotted and the data are offset for 
clarity) . The orbit assignments are indicated by symbols and 
follow the notation of Figure Q 
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clearly the increased damping with shorter mean-free- 
path I and increased average Fermi wavevector (hp) of 

the orbit. R T = A/(sinhA) where X = ^ B - ln ^ : , 
m* is the quasi-particle effective mass. Finally, the spin 
splitting factor is given by Rs = cos( 7r " 9 "^ 1+ ' 5 '^ ) where 
(1 + 5) is the orbitally averaged exchange-correlation 
(Stoner) enhancement factor, g is the electron (/-factor, 
m e is the free-electron mass and n is an integer. 

The electronic structure and dHvA orbits of MgBj 
have been calculated by three different groups [ill ll2L 
Hjj . The calculated Fermi surface is shown in Fi g. H] 
together with the expected dHvA extremal orbits |14j . 
The calculations pj show that the electronic states near 
the Fermi level arise primarily from the boron atomic or- 
bitals. The calculated Fermi surface is composed of four 
distinct sheets. Two of these arise from boron a orbitals 
and are quasi-two dimensional warped cylinders running 
along the c-direction, whereas the other two are tubular 
networks with larger c-axis dispersion, that are mainly 
formed from the boron 7r orbitals. In total 9 primary 
extremal orbits have been predicted, and 7 of these are 
labelled in the Figure (two orbits with frequencies >30kT 
have been omitted). Calculations of the dHvA frequen- 
cies and masses of the various orbits by the three groups 
are all in good agreement (the differences are typically 
100-200 T, and < 5% respectively). 

The fast Fourier transforms (FFTs) of the raw torque 
data between 20 and 32 T at 1.4 K for crystal K are 
shown in Fig. [21 as a function of angle 9 as the sample 



FIG. 3: Observed frequencies versus angle as the samples were 
rotated from H\\c to (approximately) H\\a. The solid lines are 
polynomial fits to the data as described in the text, and the 
dotted lines are guides to the eye. 



was rotated from H\\c to H\\a (H\\c = 0°). In addition 
to the frequencies observed in our previous study (i 7 ! , 
F 2 and F 3 ) several additional peaks are visible. Some 
care is needed in interpreting these new frequencies. Be- 
cause the cantilever is deflected slightly by the torque, 
the torque measurements are not made at constant an- 
gle. This generates spurious harmonics and combinations 
of the main dHvA oscillations. For crystal K only weak 
harmonics of F 2 are observed (just above F 5 - see Fig.|2J, 
whilst for the larger crystal B, several harmonics of F3 
and a frequency corresponding to F% + F2 were observed. 

The observed frequencies (omitting those assigned as 
harmonics or combinations) are shown in Fig.[3]as a func- 
tion of for both crystals. The solid lines in Fig.[3]are fits 
of the observed F(6) values to F(6) = Yli=i a<i/ cos 1 (8 — 
0o ) (^o = or 90°), which we use to extrapolate the 
observed frequencies to the symmetry points (we denote 
these frequencies as F°). The assignment of the FFT 
peaks to the orbits shown in Fig. ^ was achieved by com- 
paring the values of the frequencies obtained by extrap- 
olation, and their angular dependencies, with the calcu- 
lations. 

For crystal K, signals from 6 orbits associated with all 
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TABLE I: Summary of dHvA parameters for both samples (K and B) along with the theoretical predictions fThHl2| 
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550 0.548 ± 0.02 


0.251 


1.18 + 0.1 
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546 


380 


0.553 + 0.01 


1.20 ± 0.04 
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0.648 + 0.01 
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0.47 
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0.441+0.01 


0.40 ± 0.03 
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0.31+0.1 
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4 sheets of Fermi surface are observed. We are therefore 
able to verify experimentally the Fermi surface topology 
predicted by Kortus et al. [l| and shown in Fig. ^ F° r 
crystal B, in addition to the 3 orbits observed previously, 
3 further frequencies are seen. One of these can be as- 
signed to F A . The frequency of F 3 * (F=4600 T) is close 
to that predicted for orbit 7 (F=4294 T)HH (a nd has 
a similar effective mass), however, the frequency is also 
close to that expected for an orbit equivalent to F3 but 
from the tube oriented along a* (i.e., at 60° to a). A sub- 
sequent in-plane rotation study showed the latter to be 
the most likely origin. The origin of F^s is less clear, but 
it could arise from a slight warping of the same in-plane 
tube responsible for F3. 

Table [I] shows that the values of F° found for the 
two crystals are in good agreement. In total, we have 
studied crystals from six batches (5 from Tokyo and one 
from Zurich) and so far, the dHvA frequencies agree to 
within 30 T or 0.06 % of the basal area of the first Bril- 
louin zone (50.2 kT). This suggests that although T c is 
slightly reduced compared with the best polycrystalline 
samples, any possible Mg deficiency is very reproducible, 
and probably small [IsT ]. 

The differences between the measured F° values and 
those predicted by theory is a significant fraction of F° 
in some cases, but we note that 100 T discrepancy only 
amounts to ~ 0.3 % of the basal area of the Brillouin 
zone. The volumes of the tubes are proportional to the 
average of the two extremal areas Il3l an d these are both 
~ 16% smaller than the calculations [13, implying a cor- 
responding reduction in the number of holes in these two 
tubes 01 . This may have a significant effect on calcula- 
tions of physical properties such as T c , its dependence on 
alloying or pressure, and the London penetration depth. 

It is instructive [I2I Il3| to calculate the Fermi energy 
shift AE which would bring the theoretical frequencies 
in line with experiment. As the dHvA band mass is de- 
fined as vriQ — §^ , the necessary band shifts are given 
by AE — ^AF (where AF = F Th - F eKp ). There is 
remarkable consistency between the AE values for the or- 



bits, with the values roughly falling into two groups (Ta- 
ble [J. For the a sheet orbits (1,2,5,6) the average shift 
is 83 ± 4 meV, whereas for the it sheet orbits (3,4) it is 
—61 ± 5 meV. Because of the high degree of reproducibil- 
ity of the frequencies between samples, it is unlikely that 
this discrepancy is caused by sample impurities or non- 
stoichiometry. Instead, it seems to imply a shortcoming 
of the LDA calculations [Tlli"^ . 

Quasiparticle effective masses were determined by per- 
forming field sweeps at different temperatures and fitting 
the amplitudes to Eq. ^ Our experimental values of m* 
are compared with the calculated band masses in Table 
QJ In MgB2 the dominant source of mass enhancement is 
the electron phonon interaction. If we assume this is the 
only source of enhancement we can calculate an upper 
bound for the electron-phonon coupling constants A cp , 
from A op = m*/m B - 1. The results (Table show that 
the values of A ep on both the a sheets are approximately 
a factor three larger than those on the it sheets. 

A detailed comparison with the orbit-resolved theoreti- 
cal values ^3 °f ^ ep i s a l so shown in the table. Generally 
our values are slightly smaller than the theoretical ones. 
The most significant differences are for the larger orbits 
on the a tubes for which A op values which are ~ 20% 
smaller than theory. The small shifts in Fermi level de- 
scribed above do make small differences to the calculated 
band masses increasing the A| p p values by ~ 6%, 
but the differences remain significant. These differences 
could have a relatively large effect on the detailed calcu- 
lations of T c . For example, using the isotropic McMillan 
equation, with parameters appropriate to MgB2 we 
estimate a 20% reduction in A ep would imply a 8 K re- 
duction in T c . It is likely that the reduced A values we 
observe are caused by phonon anharmonicity, which has 
been shown to reduce the average value of A by around 
20 %0. 

Scattering rates (t _1 ) were calculated by fitting the 
raw torque versus field data to Eq. ^ The results are 
also shown in Table |U It can be seen that the mean free 
path (£) for sample K is significantly longer for orbits 
F2, F\ on the c-axis tube but not for F3 on the in-plane 
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Angle, 9 

FIG. 4: FFT amplitude versus angle for frequencies Fx, Fi 
and Fz in crystal K. The solid lines are fits to the data as 
described in the text. The dotted line is a guide to the eye. 

tubular network. The reason for this is not clear, but 
it does not necessarily indicate a higher crystal quality. 
However, F5 and F§ are almost certainly observable in 
crystal K because of the larger values of I for the two 
c-axis a tubes. These differences in I may be at least 

partially responsible for a significant difference in B C 2 

He 

values for the two crystals. For crystal B, B^ 2 — 3.3 T, 
= 17 T, and for crystal K, B^ = 2.5 T, B^f = 12 T 
(all values quoted at T=0.3 K). 

The longer mean free path of the smaller a tube (F\ 
and F2) in crystal K, has allowed us to track the ampli- 
tude of the signal over a very wide range of angle (up 
to 67° and 81° respectively). The data shown in Fig. 0] 
show pronounced minima at 9 = 51 ± 1°, 9 — 53 ±1°, 
6 = 75.7 ± 0.5° for Fi, F2, and F3 respectively, which we 
attribute to the spin-zero effect (the R s damping term in 
Eq. [JJ. Using the calculations of ms(^) by Harima 
we can deduce the enhancement of the spin susceptibil- 
ity (1 + S) for these three orbits (at the spin zero angle). 
Taking g = 2, we find that 5=0.07, 0.12 and 0.45 for the 
three orbits respectively. The reason for the significantly 
larger enhancement on the ir band is not clear. Band 
structure calculations 0] predict S — 0.31 and 0.26 on 
the F2 and F3 orbits respectively. The calculation there- 
fore overestimates S on the a sheet by a factor 2.5 and 
underestimates it on the ir sheet by a factor 1.7. 

The solid lines in Fig. 01 show a fit to the data for F\ 
and F2 with Eq.Q] using the tub values of of Harima [111 ]. 
There are 3 free parameters: the overall amplitude, r _1 
(assumed constant as a function of angle), and S. The 
fit is remarkably good and the values of I were found to 
be 410 A and 900 A for F\ and F2 respectively which 
are close to the values obtained from a fit to the field 
dependence at a single angle given in Tabled However, it 
was not possible to fit the data for F3 because the abrupt 



fall for 9 < 68° cannot be explained simply in terms of 
the angular dependence of and a constant scattering 
rate (Dingle plots show that r is approximately constant 
with angle). 

In conclusion, our data strongly support the overall 
topology of the predicted electronic structure of MgB 2 
and the calculations of the clcctron-phonon coupling con- 
stants for the different orbits. Our data give direct ev- 
idence that the electron-phonon interaction is large on 
both c-axis a sheets and much smaller on both ir sheets. 
We have therefore obtained conclusive evidence in favor 
of the two key ingredients in the two-gap model of super- 
conductivity in this compound, namely the Fermi sur- 
face topology and the disparity in the electron-phonon 
coupling for the a and n bands. 
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